Optical conductivity of strongly correlated electron systems 
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We present an exact diagonalization study of the frequency and wave vector dependent conduc- 
tivity a(q,to) in small clusters of 2D t—J model. Unlike the related dynamical density correlation 
function, o(q = 0, w) in the underdoped regime has the exchange constant J as its characteristic 
energy scale and is dominated by a resonance-like excitation with frequency ~ 1.7 J. We interpret 
this as transition to a p-like excited state of a spin-bag type quasiparticle (or, alternatively, a tightly 
bound spinon-holon pair) and show that a simple calculation based on the string picture explains 
the numerical results semiquantitatively. For doping levels > 25% t remains the only energy scale 
of o(q = 0, u>). 
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Electronic excitations at finite frequency seem to be a 
common feature of doped Mott-Hubbard insulators, and 
are observed experimentally in the optical conductivity of 
cuprate superconductors |l]] , as well as in numerical stud- 
ies of 2D strong-correlation models |2||j| . It is the purpose 
of the present paper to present a systematic study of the 
optical response function and its dependence on both, 
parameter values and hole concentration, for the stan- 
dard t—J model. This model is supposed to describe the 
low energy physics of the CuOi planes, at least for ener- 
gies which are below the 'binding energy' of a Zhang-Rice 
singlet; the latter may be estimated to be ~ leV. As a 
key result, the optical conductivity <j(q — 0, u>) in the 
underdoped regime is shown to be dominated by a single 
excitation which has the exchange constant J as its rele- 
vant energy scale, in contrast to e.g. the behaviour seen 
for the closely related (see below) dynamical charge cor- 
relation function |Q. We show that a simple calculation 
based on the string picture can explain this dominant ex- 
citation as a transition to a p-like excited state of a spin 
bag-like quasiparticle or, alternatively, a tightly bound 
spinon-holon pair. We study the finite frequency optical 
response, defined as 

a a (q,w) = V- \(^\j a (q)\^o)\ 2 S(w - (E u - E )), 



where \ip v ) (E v ) denotes the v th eigenstate (eigenenergy) 
of the system (in particluar v—Q denotes the ground 
state). Also, j a with a — x,y denotes a component of 
the current operator 



i 2J U, 



iq-{R m +R n )/2 r B 



R„ 



In the present study we restrict ourselves to the finite- 
frequency response and disregard the 'Drude peak' at 



q=Q and oj=Q. The latter is absent in finite clusters 
with periodic boundary conditions (as used in the present 
study), so that its weight D can only be inferred indi- 
rectly via a sum rule: 
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The Drude weight D and its dependence on hole doping 
and parameter values have been studied previously by 
various authors. In 2D systems and hole concentrations 
5 < 0.25 to good approximation D ~ S, with the constant 
of proportionality being nearly independent of J/t [^|. 
Since D originates from the free acceleration of charge 
carriers in an applied electric field, this result suggests 
a carrier number given by the number of doped holes, 
which is consistent with a variety of experiments nj. In 
a nonintcracting single band model the Drude part is 
the only contribution to a(q = 0,u>); by contrast, the 
finite frequency part of a(q — 0, u>) is a special feature of 
correlated systems, and the excitations which are probed 
by this part of the correlation function have not yet been 
identified. In the present study we focus exclusively on 
this feature of the frequency dependent conductivity. 
To begin with, Figure [I] shows the 'dispersion' of the 
current correlation function w ■ <j(q,uj) for different J/t. 
The 'dominant features' do not change appreciably with 
J, i.e. the energy scale of the correlation function is t (an 
exception is the largest peak in the (7r/3, 7t/3) spectrum, 
which seems to scale with J for J>1). This fact as well 
as the dispersion of the spectral weight, which resembles 
a 'smeared out' free electron band of width ~ 8t, is to be 
expected on the grounds of analogous results for the dy- 
namical charge correlation function (DCF) B. For e.g. 
q parallel to (1,0) one has [H,n q ] — 2sm(q x /2)j x (q), so 
that the current correlation for finite momentum transfer 
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q is closely related to the DCF, and in principle provides 
no new information. This is not the case, however, for 
the optical conductivity, a(q — 0,u>), so that we hence- 
forth consider only case q—0 and drop the wave vector q 
for simplicity. 
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FIG. 1. Current correlation function u-a(q, ui) for different 
momenta along the (1, 1) direction in the 18-site cluster, with 
two holes, (^-functions are replaced by Lorentzians of width 
0.5t 

Then, a(u>) is shown in Figure || for 1 and 2 holes and 
different J, for all available cluster sizes N (except for 1 
hole in 10 sites, J = 0.125, which already is in the fully 
polarized Nagaoka state) . Comparing with the spectra in 
Figure [l] which correspond to large momentum transfer, 
such as (it, 7r), some substantial differences are obvious: 
for a single hole as well as for two holes and larger values 
of J the most prominent feature is a pronounced maxi- 
mum of intensity at the lower edge of the spectra, whose 
energy scales with J. This maximum seems to have a 
more resonance-like character (i.e. a substantial 'broad- 
ening') for two holes and smaller J; even in this case, 
however, the spectra do have a J-dependent offset which 
is nearly the same as for the single hole. The 16-site 
cluster is exceptional, in that the peak still is very pro- 
nounced even for two holes, which may be an effect of 
its special geometry. As mentioned above, for the small- 
est momentum transfer, (7r/3, 7t/3), and large J a similar 
scaling with J can be seen also in Figure |l|. Since a(q, u>) 
ultimately must be a continuous function of q, there has 
to be a crossover of energy scales from small to large q. 
The wave vector where this crossover occurs apparently 
depends on J and probably is smaller than (7r/3, 7t/3) for 
J<0.5. 

Next, the evolution of the optical conductivity with dop- 
ing is shown in Figure [| There is a rather obvious 
crossover between the low doping region with a J de- 
pendent offset, and the higher doping region where t is 
the only remaining energy scale. For 'physical' values 
of J, i.e. J ss 0.4, the two-hole spectra in 16, 18 and 
20 sites(corresponding to (5=0.1 — 0.125) still show the 



J-dependent offset, wheras already at <5=0.2 (2 holes in 
10 sites) no more more scaling with J can be seen, so 
that the critical concentration S c falls between these two 
values: 0.125<i5 c <0.2. Similarly one may estimate that 
0.2<S C < 0.25 for J=l, and <5 C «0.25 for J=2. It is tempt- 
ing to speculate that the crossover for physical values is 
close to the 'optimal doping', in the the cuprates, (5=0.15 
but it should be noted that exact diagonalization allows 
to vary the hole concentration only in rather coarse steps, 
so that it is difficult to give a really accurate estimate for 
5 C . We note that such a crossover is also quite consistent 
with analogous results for spin and charge correlation 
functions j|,[7) . Moreover the doping dependence of the 
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FIG. 2. Optical conductivity a x (ui) for low hole concentra- 
tions in different cluster size N. The thin lines indicate the 
points ui=1.7J, Lorentzian broadening O.lt. 

Drude weight changes from the strong-correlation form 
D (x 8 to essentially free electron-like behaviour at ap- 
proximately the same hole concentration ||. All of this 
indicates that a profound reconstruction of the electronic 
structure takes place at this concentration. 
We focus on the physically more interesting low doping 
region and proceed to an interpretation of the dominant 
'J-peak' in (t(uj). The similarity between the single and 
two hole cases suggests that the essential physics is al- 
ready present for a single hole in an antiferromagnet, a 
limiting case that is reasonably well understood in terms 
of the 'string' picture. Following Ref. Q we assume that 
'charge spectra' in the low doping region are dominated 
by excitations of internal degrees of freedom of the spin 
bag like quasiparticles. More precisely, we assume that 
the quasiparticles in the low doping regime correspond 
to a hole oscillating rapidly (i.e. on an energy scale ~ t) 
within a region of enhanced spin disorder ||-|o). The 
entire quasiparticle (i.e. 'bare hole' plus 'defect region') 
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then drifts slowly (i.e. on an energy scale J) and coher- 
ently through the system. 
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FIG. 3. Optical conductivity a x (u}) for JV=16 and different 
hole numbers. Lorentzian broadening O.lt. 

A different way to express this picture would be 'incom- 
plete spin charge separation', where the 'holon' (=bare 
hole) and the 'spinon' (=first spin defect created by the 
hopping hole) remain tied together firmly by a 'string', 
resulting in a bound state. As discussed in Ref. Q both, 
the different energy scales of spin and charge correlation 
function and their different scaling behaviour with hole 
number rih can be understood immediately by assuming 
that the density operator couples only to the rapid inco- 
herent motion, whereas the spin operator couples only to 
the slow coherent motion. Extending this line of thought 
we would expect that the current operator, which acts on 
charge degrees of freedom as well, should also couple to 
the rapid incoherent motion, i.e. excite 'internal degrees 
of freedom' of the spin bags. Due to the high symmetry 
for q = 0, however, there is predominantly one special 
type of state which is excited by the current operator: 
we may assume, that the drifting spin bag (or, phrased 
differently, the bound spinon-holon pair) to good approx- 
imation has maximum point group symmetry, i.e. it is an 
s-like state. The current operator for q = is a vector, 
and hence must obey a 'dipole selection rule', in other 
words: it can couple only to 'p-like' excited states of the 
spin bag (or spinon-holon pair). 

We make these considerations more quantitative, using 
the formalism of Ref. || . We begin with a hole created 
at some site i in the Neel state. Acting with the hop- 
ping term, we obtain 'string states', where the hole is 
connected to its 'starting point' j by a trace of misal- 
ligned spins (ll). We denote a state with v such defects 
as |i, v, V) where V is shorthand for a set of numbers 
which parameterize the geometry of the hole path. The 



wave function for a hole which is trapped in the string 
potential around the site i then reads 

itfi,A> = x>^5>A(p)Mp>. (!) 

v V 

Here, <fi(V) denotes an extra phase factor, which deter- 
mines the 'orbital symmetry' A of the self-trapped state. 
Let e denote a unit vector in the direction of the first 
step of the path V . Then, 4>\(V) = I will give a totally 
symmetric ('s-like') spin bag, (f>\(V) — e x will give a state 
with p^-type symmetry etc. The coefficients a^\, which 
play the role of a 'radial wave function' associated with 
the orbital symmetry A, are assumed to depend only on 
the length of the string, v. If we assume that the mag- 
netic frustration in the system increases linearly with the 
string length v frH| , and introduce f3^\ — ck m ,a(z — l)^ 1 / 2 
(with z the coordination number), the latter function can 
be determined from the (numerical) solution of the 'ra- 
dial Schrodinger equation' |^ 

^-tft.A + 2 J/V = E x f3 ,x 

z — 1 

-t( A,+1,A + &-1.A ) + J{\ + v)p VlX - £ A A,,A, (2) 

where t = y/z — It. The normalization condition reads 

with n s = 1, n p = 2 and /3o,p = 0. (and consequently 
the first equation (^) being omitted for A = p). Next, a 
propagating spin bag with orbital symmetry A and mo- 
mentum ko would be described by 

When comparing with the numerical results we choose 
the 18-site cluster, where the single hole ground state 
momentum is fco=(27r/3, 0) (this avoids complications 
due to the spurious degeneracies in the 4x4 cluster). 
The little group of feo comprises E and the reflection by 
the (1, 0) direction, the ground state wave function being 
even. Choosing the current in (1,0) direction the final 
states must be even as well and thus should correspond 
to a propagating p^-like spin bag, whereas choosing the 
current in (0, 1) should couple to a propagating pj,-like 
spin bag. Assuming that the dispersion of the states (||) 
is solely due to the relaxation of the strings through the 
transverse part of the Heisenberg exchange and follow- 
ing Ref. H in the computation of the respective matrix 
elements, we obtain the dispersion of the s-like state as 
E s (k) = E s + 4h s ([cos(k x ) + cos(k y )] 2 — 1), and that 
of a p a state as E p (k) = E p + 2h p cos(2k a ), with the 
'spin-flip matrix element' h x = S^Lo /^aA,+2.a- 

The latter describes the truncation of the string by the 
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Heisenberg exchange. Differences between those ener- 
gies at feo give the excitation energies in the optical 
conductivity. The matrix element of the current op- 
erator is g(fco)b'a|* Pa {ko))=-2itY,7=o( a V;S a v+i,p ~ 
otu,pOt v +i,s)- Then, Figure ^ compares the current corre- 
lation function obtained by Lanczos diagonalization with 
the result of the string calculation for transitions to the 
lowest p-like state. 
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FIG. 4. Exact current correlation function for one hole in 
18 sites (fulll line) compared to the result of the string calcu- 
lation (dotted line). Lorentzian broadening 0.025t. 

We note that the string calculation correctly repro- 
duces the J-dependence of the offset. For larger J and 
current perpendicular to the momentum, even the weight 
of the dominant peak is reproduced quite well. On the 
other hand, for smaller J the 'peak' at the bottom of 
the Lanczos spectra seems to 'decay' into a multitude of 
small peaks, i.e. there is a strong damping. This is to 
be expected, because the p-like state first of all is high 
in energy (on the energy scale of spin fluctuations) and, 
having a node in the center, necessarily has a large spa- 
tial extent; most probably it is therefore very susceptible 
to scattering from spin excitations and/or finite size ef- 
fects, so that the p-like 'state' rather has the character 
of a 'resonance'. While it seems that the string calcula- 
tion reproduces the rapid decrease of spectral weight, any 
more quantitative comparison must fail due to this decay 
of the peak. Taking into account its highly oversimplified 
nature, the string calculation nevertheless predicts some 
trends reasonably well, and we believe that it actually 
is the correct explanation for the observed behaviour. It 
should be noted that a calculation for 'spinon'-'holon' 
pair, which is bound by a potential ~ J times their dis- 
tance would almost inevitably give an identical result. 
Finally, the apparent scaling of the dominant peak in 
the current-correlation function for q=(jr/3, 7r/3) with J, 
which seems to occur for larger J, may be an indication 
that a propagating p-like spin bag represents an approx- 
imate eigenstate also for small finite momentua. 
In summary, we have presented a systematical study of 
the optical conductivity in the 2D t— J model. We found 



that over the entire 'underdoped' range of hole concen- 
trations, the optical spectrum is dominated by a single 
excitation at an energy of w 1.7 J. The appearance of the 
exchange constant J as the energy scale of a 'charge spec- 
trum' is not consistent with an interpretation in terms of 
spin charge separation, where the charge spectra would 
be described by the excitations of holons, whose energy 
scale is t; the latter would be neccessary to be consis- 
tent with the finite momentum current correlation and 
density correlation functions. On the other hand, an in- 
terpretation in terms of the 'string' picture, (or alterna- 
tively the assumption of tightly bound spinon-holon pair) 
where the dominant excitation in tr(w) is interpreted as 
transition to the lowest excited state of the 'spin bag' 
with p-type symmetry, readily provides a semiquantita- 
tive explanation for the numerical data. The absence 
of an appreciable finite frequency component in ID then 
fits nicely with this picture: in ID spinon and holon truly 
unbind, so that excited 'bound states' which could serve 
as final states in <j(uj) do not exist. 

We note that 'incomplete spin-charge separation', i.e. the 
formation of tightly bound spinon-holon pairs, in reality 
would imply absence of spin charge separation. When 
probed with wave lengths larger than its spatial extent, 
or very low energies, the spinon-holon pair should be- 
have just like a spin-1/2 Fermion, so that the 'effective' 
low energy theory would be a Fermi liquid of quasipar- 
ticles corresponding to the doped holes. This scenario 
is indeed suggested by diagonalization [jl2],[l3| and other 
p4] , p^[ studies. The nature of the 'transition' at higher 
doping levels, and the resulting ground state remains to 
be clarified. 
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